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چکیده 


اقلیم یکی از اساسی‌ترین عوامل طبیعی است که تغییر آن از مهم‌ترین چالش‌های قرن جاری می‌باشد. انتشار روز افزون گازهای گلخانه‌ای به عنوان 
مهم‌ترین عامل تأثیرگذار بر تغییر اقلیم منجر به تغییر دماء بارندگی و دیگر پارامترهای اقلیمی می‌شود. بخش کشاورزی در مقایسه با سایر بخش‌ها بیشتر 
در معرض آسیب‌های ناشی از تغییر اقلیم است به طوری که تغییر الگوی بارش و دمای متوسط جو بر تولید انواع محصولات باغی و زراعی که 
عمده‌ترین منابع غذایی کشور را تشکیل می‌دهند. آسیب وارد می‌کند. مدیریت ریسک در بخش کشاورزی به ویژه ریسک آب و هوا از اهمیت خاصی 
برخوردار است اما استراتژی مدیریت ریسک بدون شناسایی و اندازه‌گیری آن نمی‌تواند به درستی اجرا شود. از این رو مطالعه حاضر درصدد است ریسک 
عملکرد ناشی از تغییر اقلیم گندم دیم شهرستان اهر در دو دوره پایه (۱۳۹۴- ۱۳۶۵) و آتی (۱۳۹۵-۱۴۲۴) را با بکارگیری مدل ارزش در معرض خطر 
(۷۵8) و ارزش در معرض خطر شرطی (۷۸18)) آب و هواء مدل آکواکراپ مدل گردش عمومی جو 11۸0613 تحت سناریو انتشار ۸2 و ریزمقیاس 
نمایی 175-2 اندازه گیری کند. برای این منظور اطلاعات هواشناسی در دوره (۱۹۸۶-۲۰۱۵) و عملکرد سال‌های زراعی ۱۳۷۰-۹۴ شهرستان اهر 
گردآوری شد. نتایج مطالعه نشان می‌دهد که در دوره آتی متغیرهای بارش (به میزان ۹/۱۵ میلی‌متر)» حداکثر دما (۰/۶۳ درجه سانتی‌گراد4 حداقل دما 
(۰/۵۶ درجه سانتی‌گراد» تابش خورشیدی (۰/۰۶ مگاژول بر متر مربع در روز)» حداکثر رطوبت نسبی (۰/۶۱- درصد) و حداقل رطوبت نسبی (۰/۲۹- 
درصد) سرعت باد (۰/۰۳ متر بر ثانیه) و تبخیر و تعرق (۱۸/۴۸ میلی‌متر بر سال) به مقادیر فوق تغییر پیدا می‌کنند. براساس تغییر متغیرهای آب و هوایی 
متوسط عملکرد در دوره پایه از ۰/۹۵۴ به میزان ۰/۹۹۹ تن در هکتار در دوره آتی افزایش یافت. برای اندازه‌گیری ریسک توزبع بتا و ویکبای برای 
عملکرد مناسبترین توزیع انتخاب شد. مقدار کاهش ریسک عملکرد در سطوح احتمالاتی ۱ ۵ و ۱۰ درصد در روش ۷۵ برای توزیع بتا ۰/۱۳۹ ۰/۸۳ و 
۱ برای توزیع ویکبای ۰/۲۹۷ ۰/۸۱ و ۰/۳۱ تن در هکتار و در روش 0۷۵8 در توزیع ۰/۱۵۸ ۰/۱۱۵ و ۰1۹۳ و برای توزیع ویکبای ۰/۴۰۳ و 
۸ ۷۷ من :در هکتار پزست امد پراسانن نتایم,مطالبه تبخیرو فرق یکی از مشیرهانی ات که موجب آفزایفن ویشکه عملکرد می شود از ای 
رو برای کاهش تبخیر و تعرق از سطح خاک استفاده از خاک‌پوش‌های طبیعی و مصنوعی به کشاورزان توصیه می‌شود. 


واژه‌های کلیدی: ارزش در معرض خطر» اهر تغییرات اقلیم» ریسک آب و هو گندم دیم 


مقدمه 


تولید در بخش کشاورزی با سایر بخش‌های اقتصادی تفاوت‌هایی 


۰ ۲ و ۲- به‌ترتیب استاد» کارشناس ارشد و دانشیار گروه اقتصاد کشاورزی» 
دانشکده کشاورزی» دانشگاه تبریز 
(*+- نویسنده مسئول: .6۵۱۷۵۵0 1)طافهلداع :اتفعصط) 
۴- دانشیار گروه مهندسی آب دانشکده کشاورزی, دانشگاه تبریز 
102025 :]10 


قرار می‌دهد عوامل طبیعی نیز تولید و درآمد در اين بخش را به شدت 
با چالش مواجه می‌سازد. اقلیم را می‌توان از اساسی‌ترین عوامل 
طبیعی دانست که تغییر آن یکی از مهم‌ترین چالش‌های قرن جاری 
است. انتشار روزافزون گازهای گلخانه‌ای به عنوان مهم‌ترین عامل 
تأثیرگذار بر تغییر اقلیم منجر به تغییر دماه بارندگی و دیگر متفیرهای 
اقلیمی شده و در نهایت باعث تغییر در میزان رطوبت خاک افزایش 
سطح آب 

درياهاء افزایش دما و وقوع سیل و خشکسالی در بسیاری از 
مناطق دنیا شده است (۱۲). از این رو انتظار می‌رود بخش کشاورزی 
به عنوان تأثیرپذیرترین بخش اقتصادی به تغییرات آب و هوا در 
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معرفن این تخطرات قرار بگیره: (۳۰) تالیرانی که:ریسک ناشی از 
تغییر اقلیم انتظار می‌رود در این بخش داشته باشد شامل تغییر در 
عملکرد و کیفیت محصول به دلیل دمای بالا و تغییر در الگوی 
بارندگی است (۲۴). انتظار می‌رود در آینده نیز افزايش دما تغییر در 
مقدار و الگوهای بارش و افزايش در تکرار و شدت مخاطرات محیطی 
مات کش کسال ۶ گورههای تایبا ون شمه ۳ 

هر گونه تغییر در مولفه‌های بارش و دماء بر میزان تبخیر- تعرق 
نیز تأثیر خواهد گذاشت در نتیجه کیفیت و عملکرد محصولات 
کشاورزی در معرض ریسک تغییر اقلیم قرار خواهند گرفت. همچنین 
به دلایل ریسکی بودن فعالیت‌های کشاورزی و کاهش عملکرد ناشی 
از تغییر اقلیم انتظار می‌رود سرمایه‌گذاری و صادرات در این بخش 
کاهش و واردات و وابستگی کشور به سایر کشورها افزایش یابد که به 
نوبه‌ی خود. تأثیر منفی بر تراز تجاری اقتصاد دارند. از اين رو معمولا 
برای اجتناب و يا به حداقل رساندن آثار ریسک آب و هوایی در بخش 
کشاورزی اقداماتی انجام می‌گیرد که می‌توان به صورت زیر 
گروه‌بندی کرد: ۱- اقدامات مدیریتی ۲- اقدامات فنی و تجهیزات ۳- 
اقدامات سیاستی و اقتصادی ۴- اقدامات زیر بنایی (۲۴). اقدامات 
سازگاری که در بالا ذکر شد تا حدودی عمومی است و که برای درک 
تأثیرات اقدامات نیاز به کمیت ریسک تولیدات کشاورزی دارد. ولی به 
طور کلی ریسک آب و هوایی ناشی از تغییر اقلیم در تقسیم‌بندی 
خطرات بخش کشاورزی در قسمت شدت پایین خسارت با فروانی بالا 
قرار می‌گیرد (۲۳)» همچنین این نوع ریسک جزو ریسک‌های غیر 
مصیبت‌بار آب و هوا کمتر مورد توجه است به طوری که زیان‌های 
ناشی از ریسک‌های مصیبت‌بار توسط بیمه محصولات کشاورزی 
پوشش داده می‌شود اما برای ریسک‌های غیر مصیبت‌بار فقط اقدامات 
غیر مصیبت‌بار را دو چندان می‌کند. با اين وجود به دلیل وابستگی 
تولیدات کشاورزی به آب و هوا ریسک تغییر اقلیم تأثیر بسزایی در 
این ریسک نامحسوس کمی شود تا در اقدامات مدیریتی بیش از پیش 
مورد توجه قرار گیرد. 

براساس مطالعات صورت گرفته در داخل کشور دمای کشور تا 
محصولات کشاورزی نیز در معرض ریسک خواهند بود که براساس 
کرش ارم اطع معط ابیت رم قوف فرلند لاش 
کشور (که با سطح زیر کشت معادل ۷۳/۱ درصد عمده‌ترین نوع 
مقر لاک ماکان کشور میاه تاسال ۲۰۲۷ صقر ای ۸ 


رارسا ی ۵ ای ها هشال ماب ۸ 
الی ۳۰ درصد کاهش یابد (۷). 

در بین گروه غلات نیز گندم در سال زراعی ۱۳۹۲-۹۳ با 
بیشترین سطح زیر کشت معادل ۵۱/۲۰ درصد مهمترین محصول 
کشاورزی کشور بوده است. استان آذربایجان شرقی با سطح زیر کشت 
۰ درصد گندم کل کشور رتبه ششم در تولید این محصول را به 
خود اختصاص داده است که ۷۸/۷ درصد سطح زیر کشت استان به 
گندم دیم و ۲۱/۳۲ درصد آن به گندم آبی اختصاص دارد (۲۲). در 
میان شهرستان‌های استان, اهر به عنوان یکی از بزرگترین 
تولیدکنندگان گندم» ۵/۹۵ درصد سطح زیر کشت استان را بخود 
اختصاص داده که در این بین» سهم اراضی دیم ۸۹/۲۵ و اراضی ان 
۵ درصد می‌باشد (۱۵). تغییرات بارندگی و دما ریسک‌های آب و 
هوایی را برای گندم دیم نشان می‌دهند که مهم‌ترین متفیرهای 
نشان‌دهنده تغییر اقلیم هستند. معمولاً بیشترین میزان وابستگی بین 
عملکرد گندم و میزان بارندگی اتفاق می‌افتد به گونه‌ای که ضریب 
همبستگی بین عملکرد گندم شهرستان اهر و بارندگی در فصل رشد 
گندم (مهر تا تير ماه) ۰/۴۵ بدست آمد که نشان می‌دهد بیشتر 
تغییرات عملکرد توسط بارش توضیح داده می‌شود. از این رو انتظار 
می‌رود پدیده تغییر اقلیم بر تولید گندم دیم در شهرستان اهر تأثیر 
به‌سزایی داشته باشد که هدف مطالعه حاضر در این راستا می‌باشد. 

تاکنون مطالعات کسترده‌ای در خضومن بررسی تالیر تغییر انیم 
بر تولید محصولات کشاورزی در داخل و خارج صورت گرفته است. 
کوچکی و کمالی (۱۹) در تحقیقی با ستفده از مدل گردش عمومی 
(00) میزان کاهش گندم دیم در اثر تغییر اقلیم را در کشور بین 
دو دوره ۲۰۲۵ و ۲۰۵۰ شبیه‌سازی کردند. نتایج شبیه‌سازی نشان داد 
که شدت کاهش عملکرد گندم دیم در مناطق دیم‌خیز کتسور برای 
سال ۲۰۲۵ در محدوده ۱۶ تا ۲۴ درصد و برای سال ۲۰۵۰ میلادی 
در دامنه ۲۲ تا ۳۲ درصد می‌باشد. در مطالعه‌ای دیگر آبابایی و 
همکاران (۱) تغییرات اقلیم و تأثیر ریسک ناشی از آن را بر عملکرد 
گندم با استفاده از مدل ٩۹۲۷۸۵۳‏ در منطقه رودشت اصفهان برای 
دوره ۲۰۱۱-۲۰۳۰ مورد ارزیابی قرار دادند. نتایج نشان داد که متوسط 
عملکرد سبی محصول قحت دو سناریو تغییر آقلیم ۸2 و 1ظ نسبت 
به سناریو مبنا (۱۹۶۱-۱۹۹۰)» به ترتیب ۱/۴۹ و ۲/۱ درصد و متوسط 
عملکرد دانه گندم به ترتیب ۳/۱۹ و ۱۷/۹ درصد کاهش خواهند 
یافت. همچنین بسکابادی و همکاران (۵) اثر یی اقلیم ناشی از دو 
برابر شدن غلظت ۲۶ را که بر دما و بارندگی تأثیر می‌گذارده بر تولید 
گندم در دوره زمانی ۱۳۷۹-۱۴۳۰۰ با در نظر گرفتن عوامل اقتصادی 
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اندازه‌گیری ریسک عملکرد ناشی از تغییر اقلیم در گندم دیم شهرستان اهر ... ۱۴۱ 


با استفاده از مدل برنامه‌ریزی غیرخطی در مشهد مورد بررسی قرار 
دادند. براساس یافته‌های پژوهش عملکرد بهینه گندم از ۱۳۳۱/۴۵ 
سال هدف تغییر کرده است. همچنین پیش‌بینی می‌شود ۲۲ درصد 
مصرف آب و ۸۱ درصد مصرف کود در سال ۱۴۳۰۰ نسبت به سال پایه 
بیشتر شود. پیش بهار و همکاران (۲۷) آثار تغییرات اقلیمی را بر 
عملکرد ذرت دانه‌ای در ایران با استفاده از مدل ریکاردین و روش 
اقتصادسنجی فضایی در سه اقلیم گرم معتدل و سرد طی دوره 4۱- 
۰ بررسی نموده و نتیجه گرفتند که شدت نوسانات تغییرات 
اقلیمی در هر سه اقلیم مطالعه شده به اندازه‌ای بوده است که به 
عنوان عامل‌های ریسک سیستماتیک شناسایی شوند. از بین مطالعات 
خارجی بلانس (۳) تأثیر تغییرات اقلیم را در عملکرد چهار محصول 
مهم (ارزن» ذرت» سورگوم و منهوت) در ۲۷ کشور جنوب صحرای 
آفریقا طی دوره زمانی سال‌های ۱۹۶۱-۲۰۰۲ ضمن استفاده از یک 
مدل گردش عمومی جو (0)۷15)) ارزیابی کرد. نتایج مطالعه مبنی بر 
تغییر عملکرد تا سال ۲۱۰۰ در محدوده -۷۱٩‏ درصد تا ۶ درصد برای 
برای سورگوم می‌باشد و عملکرد منهوت بدون تغییر مانده است. 
همچنین پریتین سلیر و همکاران (۲۰) ریسک اقتصادی حاصل از 
نوسانات آب 9 هوا ر در بین دو دوره پایه (۰ ۰ 0۳۹۰ 9 آینده 
(۲۰۳۱-۲۰۷۰) در دو بخش صنعت گردشگری و کشاورزی (تولید 
گندم) ساردینیا ایتالا با روش ارزش در معرض خطر آب و هوا محاسبه 
و مقایسه کردند. نتایج مطالعه حاکی از آن است که ریسک درآمدی 
بخش کشاورزی در هر دو سناریو دوره پایه و آینده بیشتر از بخش 
گردشگری است: نیلبرد و همکاران (۲۴) نیز تثیرات تغیر قلیم را در 
مورد گندم دیم جزیره ساردینیا یلا و گوجه فرنگی تونس با استفاده 
از دو مدل ارزش در معرض خطر و آکواکراپ" در بافت‌های مختلف 
خاک منطقه در دو دوره پایه (۱۹۷۰-۲۰۱۰) و آینده (۲۰۴۰-۲۰۷۰) 
مدل‌سازی کردند. براساس نتایج حاصله در بیشتر بافت‌های خاک لوم 
شنی و لوم رسی تولید هر دو محصول به طور قابل توجهی در شرایط 
اقلیم آینده کاهش می‌یابد در صورتی که تغییرات در مدیریت 
کشاورزی از قبیل تغییر در تاریخ کاشت گندم همزمان با تفییر 
الگوهای بارندگی یا مالچ پاشی گوجه فرنگی در کاهش تلفات 

مرور مطالعات فوق بیانگر آن است که تفییر اقلیم یکی از 
مهم‌ترین چالش‌های قرن جاری بوده و باعث تغییر پارامترهای آب و 
هوایی از جمله دما و بارش می‌شود که تغییر هر یک از این پارامترها 
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تأثیر بسزایی در عملکرد محصولات کشاورزی از جمله محصولات 
دیم دارد. با توجه به اينکه گندم مهم‌ترین محصول کشاورزی است 
مدیریت ربسک آب و هوا در اين مورد از اهمیت بالایی برخوردار است 
این در صورتی است که مدیریت ریسک بدون محاسبه دقیق آن 
نمی تواند به درستی اجرا شود از این رو در این مطالعه سعی شده از 
یکی از سنجه‌های متداول ریسک برای محاسبه دقیق ریسک عملکرد 
ناشی از تغییر اقلیم با استفاده از مدل ۳۳1206۷13 تحت سناریو انتشار 
2 در دو دوره پایه (۱۳۹۴- ۱۳۶۵) و آتی (۱۳۹۵-۱۴۲۴) گندم دیم 
شهرستان اهر استفاده شود که می‌تواند در برنامه‌ریزی‌ها و طراحی 
سیستم‌های مدیریتی بسیار موثر باشد. بدین ترتیب هدف مطالعه 
اف آنلاژه یر شک مهد اف ار راکفا دب 
شهرستان اهر می‌باشد. 


مواد و روش‌ها 

روشی که در اين مطالعه برای محاسبه ریسک عملکرد ناشی از 
تغییر اقلیم در نظر گرفته شده است به این صورت است که ابتدا روند 
تفر اتتتفیرهای ابو فان با انشفاده اززیک ازیو اقلنمی سنا 
3 تحت سناریو انتشار ۸2) و مدل ریزمقیاس نمایی -1275 
8 بررسی و سپس عملکرد محصول با استفاده از مدل آکواکراپ 
شبیه‌سازی و در نهایت با استفاده از مدل آرزش در معرض خطر و 
ارزش در معرض خطر شرطی آب و هوا ریسک عملکرد محاسبه 
می‌گردد. از این رو اين بخش در سه قسمت الف) تولید سناریو اقلیمی 
ب) شبیغسازی غملکرد وچ) نازهگیری ریسک ارآنهقده است, 

الف) تولید سناریوهای اقلیمی: محققان معمولاً برای بررسی تغییر 
یهار سار وهای اقلیس اتانوس کنیا استهاده ار 
سناریوهای انتشار و مدل‌های اقلیمی شکل می‌گیرد. سناریوهای 
اقلیمی برای پیش بینی وضعیت اقلیم در آینده از داده‌های سناریوهای 
انتشار استفاده می‌کنند و ابزار این پیش‌بینی مدل‌های اقلیمی هستند 
(۱۳). ستاریو انتشار حاویی اطلاعاتی از وضعیت اقتصادی-اجتماعی و 
میزان انتشار گازهای گلخانه‌ای در اتمسفر کره زمین است. در این 
تاویوها مقادیر کاوهای خلخاندلی با ترضی بت فاسال ۲۷۰ 
افزایش می‌یابد. هر کدام از این سناریوها مربوط به یکی از 
خرهاش تور مرن عورش باق ادن سر ام 
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3 ۱۷۵۹108 
۳- سناریو انتشار ۸۰2 در برگیرنده‌ی توسعه اقتصادی منطقه محور رشد جمعیت به 
صورت پیوسته و تقویت نیروهای جمعیتی منطقه‌ای با تأکید بر ارزش خانواده و 
رسوم خانوادگی می‌باشد. 


۱۳۲ نشربه اقتصاد و توسعه کشاورزی. جلد ۰.۳۲۲ شماره ۲. تابستان ۱۳۹۷ 


اقلیمی مدل‌های گردش عمومی جو (00۷1) به نمایندگی از 
فرآیندهای فیزیکی در جو اقیانوس سطح یخبندان و زمین معتبرترین 
و پیشرفته‌ترین ابزا, برای شبیه‌سازی پاسخ سیستم جهانی آب و هوا 
به افزایش غلظت گازهای گلخانه‌ای می‌باشد. مدل مد نظر در این 
مطالعه 11800763 تخت سناریو انتشار ۸2 می‌باشد. مدل 
3 ار مرکز تحقیقاتی هدلی برای پیش‌بینی و انجام 
پژوهش‌های آب و هوا 16۳18" در نگلیس اجرا می‌شود که به 
دلیل توقای سا بالاهر آخلب:مطااات داخی و ترجه کتار 
گرفته شده و دارای مزایایی از جمله قدرت تفکیک بالای مولفه‌ی 
آقیانوسی» هماهنگی خوب بین مولفه‌های جوی و اقیانوسی آن و عدم 
نیاز به تنظیمات شار سطحی برای بهبود شبیه‌سازی است (۱۷). 
سناریو انتشار ۸2 دربرگیرنده‌ی توسعه اقتصادی منطقه محور رشد 
جمعیت به صورت پیوسته و تقوبت نیروهای جمعیتی منطقه‌ای با 
تأکید بر ارزش خانواده و رسوم خانوادگی می‌باشد (۳۸). 

فقس غیده ممل‌های فلس فر قیرت کیک مکانی کنو انیت 
است از این رو این مدل‌ها قابل استفاده در مطالعات کاربردی با ابعاد 
کوچکتر از خود نیستند. برای مثال اغلب مطالعات هیدرولوژیکی با 
فرآیندهای کوچک مقیاس سر و کار دارند که مقیاس آنها بسیار 
کوچکتر از مقیاسی است که مدل‌های گردش عمومی جو به ما 
می‌دهند. برای حل این مشکل از روش‌های ریز مقیاس نمایی استفاده 
می‌شود (۱۱). روش ریز مقیاس نمایی مد نظر در این مطالعه مدل 
1279-۷2 می‌باشد که از روش‌های ریزمقیاس نمایی آماری می‌باشد 
که جهت شبیه‌سازی حداقل و حداکثر دما ([6*)» بارش (17777) و 
تشعشعات خورشیدی ( ۸۸/1۴09 ) در یک مکان تحت شرایط 
اقلیم حال و آینده استفاده می‌شود. 

مدل 1275-172 براساس سری‌های مولد آب و هواست که از 
توزیعات نیمه تجربی (86۳06۳0۷) برای سری‌های طول روز خشک 
و مرطوب. بارش روزانه و تابش خورشیدی روزانه استفاده می‌کند. 
توزیع نیمه تجربی یک هیستوگرام با ده بازه به صورت رابطه (۱) 


است: 
(( . (10می2یلی 00‏ ورم,ره,نه) - 7۸۳ 
)۲ 6 ک رز ۳3 


(,۵,,.م] نشان‌دهنده بازه‌ها است و ,و نشان دهنده تعداد 
رخداد مورد نظر از داده‌های مشاهده شده در فاصله نام است. توزیع 


تجربی برای هر متغیر به ۲۱ پارامتر نیاز دارد ۱۱ پارامتر در مرزهای 
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ارزش‌های تصادفی از توزیعات نیمه تجربی با انتخاب نخست یکی از 
فاصله‌ها و پس از آن انتخاب یک ارزش در بازه از توزیع یکنواخت 
بدست می‌آید و ازآن جایی که توزیع نیمه تجربی مورد نظر برای هر 
متغیر انعطاف‌پذیر است. می‌تواند طیف گسترده‌ای از اشکال را با 
تنظیم فواصل تقریب کند. این بازه‌ها براساس خواص مورد انتظار 
متغیرهای آب و هوا انتخاب می‌شوند (۳۷). 

تولید داده از طریق مدل 1275-7۷2 در سه مرحله صورت می‌گیرد 
که شامل: کالیبره کردن داده‌هاء ارزیابی 9 تولید داده‌های هواشناسی 
برای دوره‌ی آنی می‌باشد. در مرحله کالیبره کردن مدل به فایلی نیاز 
شامل داده‌های روزانه دمای حداقل و حداکثر, بارش و تابش 
خورشیدی است که مدل براساس آنها برای دوره پایه اجرا می‌شود. در 
مرحله ارزیابی» داده‌های تولید شده مدل و داده‌های واقعی (مشاهده 
شده) موجود در دوره پایه با استفاده از روابط آماری ارزیابی می‌گردند 
به طوری که برای ارزیابی اختلاف میانگین ماهانه بین داده‌های 
واقعی و داده‌های تولید شده از «آزمون ۳7» توزیع داده‌های تولید 
شده و توزیع واقعی متفیرها از آزمون «کولم وگروف- اسمیرنف دو 
نمونه‌ای « ۱ 5 -1) و برای بررسی داده‌های تولید شده با داده‌های 
واقعی از آماره‌های «مجذور میانگین مربعات خطای استاندارد۵» 
(۰)/۸/57 «مجذور میانگین مربعات خطای استاندارد نرمال شده » 
(۷۸/۹۳) و «ضریب همبستگی پیرسون » ( 7 ) استفاده می‌شود. 
پس از قسمت ارزیابی در صورت اطمینان از نتایج اقدام به شبیه‌سازی 
داده‌های اقلیمی برای دوره آنی می‌گردد. نیاز اساسی مدل در این 
مرحله داده‌های دوره پایه و خروجی یکی از مدل‌های گردش عمومی 
جو می‌باشد. در حقیقت مدل با استفاده از رفتار اقلیم در دوره پایه و 
ریز مقیاس نمایی آماری داده‌های یک مدل گردش عمومی جو 
پارامترهای اقلیمی آینده را در حد مقیاس منطقه مورد مطالعه قرار 

ب) شبیه‌سازی عملکرد: پس از شبیه‌سازی داده‌های آب و 
هوایی نیاز به مدلی داریم که با استفاده از اطلاعات آب و هوایی متغیر 
عملکرد را شییه‌سازی نماید که برای این منظور از سل آکواکراپ 
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۷- مدل 1215-1۷86 برای ایجاد داده‌های بارش» حداقل و حداکتر دما و تابش 

خورشیدی استفاده می‌شود در صورتی که در فرآیند شبیه‌سازی عملکرد نیاز به 

محاسبه تبخیر- تعرق مرجع وجود دارد که علاوه بر متفیرهای ذکر شده به 

متفیرهای سرعت باد و رطوبت نسبی حداقل و حداکثر نیز نیاز است که برای ایجاد 
این متغیرها از مدل 21117061610) استفاده می‌شود. 


اندازه‌گیری ریسک عملکرد ناشی از تغییر اقلیم در گندم دیم شهرستان اهر ... ۱۴۳ 


استفاده می‌شود. مدل آکواکراپ یکی از کاراترین و جدیدترین 
مدل‌های گیاهی است که توسط متخصصان برجسته نهادها و سازمان 
های دولتی از سراسر جهان از طرف سازمان جهانی کشاورزی و غذا 
مطرح شد. این مدل با توجه به مقدار آب مصرف شده توسط گیاه 
میزان عملکرد محصول را با یک گام روزانه برآورد می‌کند رابطه آن 
به صورت زیر است: 

(۳ (,/۳ < ۲۷ ,وک ۷1 ۲ 


که در آن ۷ عملکرد محصولء 77 شاخص برداشت رگ 
چقدر به یک نزدیکتر باشد اثر محدودکنندگی دما کمتر است» "۲۷ 
پارامتر بهره‌وری آب (ماده خشک قسمت‌های هوایی گیاه ( *5 ) را 

2 ۱ 

بیان می‌کند که در واحد سطح زمین ([ 17) از هر واحد آب تعرق 
یافته (722722) تولید می‌شود که برای آب و هوا و گیاه معین ثابت در 
نظر گرفته می‌شود. پارامتر بپهره‌وری آب از طریق غلظت دی اکسید 
کربن موجود در جو با ضرب کردن یک فاکتور تصحیح نرمال می‌شود 
به این صورت که اگر غلظت دی اکسید کربن بالای ۲۶۹/۴۱ 00۲0 
باشد فاکتور تصحیح بزرگتر از یک است و برای غلظت‌های پایین‌تر از 
آن کوچکتر از یک است) و 717 تعرق محصول می‌باشد که طبق 
رابطه (۴) بدست می‌آید: 
(۴) م2( 60 ) رمک > 7۳ 

,1669 ضریب تعرق گیاهی" متناسب با پوشش کانوپی سبز ( 
"6 می‌باشد که براساس پوشش گیاهی شبیه‌سازی شده در طول 
رفتل گیاه کنظیم می‌شوده رف ضزیت نفزم آب. (مقدار آن نین ضسفز 
و یک متفیر است که هر چقدر به یک نزدیکتر باشد اثر 
محدودکنندگی آب کمتر است) م7 تبخیر و تصرق مرجع" که در 
مدل آکواکراپ از معادله فائو پنمن- مانتیث " بدست می‌آید" که اکثر 
متفیرهای آب و هوایی از طریق رابطه تبخیر و تعرق در عملکرد گیاه 
موثر هستند (۳۱). 

ج) اندازه‌گیری ریسک: پس از شبیه‌سازی عملکرد نیاز به مدلی 
داریم که بتوان ریسک عملکرد را اندازه‌گیری نماید. روش ارزش در 
معرض خطر یکی از مهم‌ترین سنجه‌های ریسک نامطلوب است که 
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در زمینه‌های مختلف مخصوصاً در موسسات مالی و صنعت بیمه به 
کار گرفته می‌شود امروزه با توجه به گستردگی و پیچیدگی کسب و 
نها زک هایی تیامی گازها را یت کی هرت ی اس هی 
ریسک‌ها و مدیریت آنها می‌تواند تصمیم‌گیری را بهبود بخشیده و 
ارزش بنگاه‌ها را حفظ نماید. در کسب و کارهای کشاورزی یکی از 
مهم‌ترین ریسک‌هایی که بنگاه‌های تولیدی را تهدید می‌کند تغییرات 
آب و هوا است که در این زمینه نیز می‌توان رهیافت ارزش در معرض 
خطر آب و هوا را به کار برد. ارزش درمعرض خطر آب و هوا (-۷) 
۴ تنشان‌دهنده حداکثر خسارت مورد انتظار ناشی از شرایط نامساعد 
آب و هوایی در یک دوره زمانی و برای یک فاصله اطمینان معین 
است (۳۰). در علم آمان ارززش در معرض خطر (۷۵1) چارک یک 
طرفه توزیع مشاهده شده در یک افق زمانی خاص می‌باشد» به طور 
کلی برای برآورد ارزش در معرض خطر به توزیع متغیر تصادفی مورد 
نظر و دو پارامتر سطح اطمینان و افق زمانی نیاز است. که در آن 
ارزش در معرض خطر میزان ریسک را در دوره مشخصی از آینده 
اندازه‌گیری می‌کند به گونه‌ای که می‌توان براساس آن سیاست‌های 
مناسبی را اتخاذ کرد. در واقع موضوع مهم این است که باید توزیع 
احتمال تجمعی مناسب را برای متغیر تصادفی مورد نظر بدست آورد 
که ممکن است توزیع نرمال نداشته و از توزیعات «دنباله‌ی پهن»" و 
«کشیدگی بالا» تبعیت کند. اگر  )0,1(‏ 7 باشد ارزش در 
معرض خطر در سطح احتمال 0 به صورت زیر تعریف می‌شود: 
(۵) 0 > (۲۷۵۵> 2۲ -(۳)/)2۲ 
(۶) (۵) ,۳ -(5)26 < ۷۵ 

در رابطه بالا "17 وارون تابع 17 و تابع کوانتیل و (5026 ارزش 
متغیر تصادفی در میانگین شرایط آب و هوایی است (۳۹ و ۲۰). 

ارززش در معرض خطر روش‌های گوناگونی برای اندازه‌گیری دارد 
که می‌توان در سه دسته زیر تقسیم‌بندی کرد: الف) روش‌های 
پارامتری که فرض خاصی را بر توزیع متغیر مورد نظر تحمیل 
می‌کنند. ب) روش‌های ناپارامتری هیچ فرضی را بر توزیع متغیر مورد 
نظر تحمیل نمی‌کنند و از آخرین توزیع تجربی متغیر بهره می‌گیرند و 
بر این فرض استوارند که آینده نزدیک تا اندازه زیادی به گذشته 
نزدیک شبیه است. ج) روش‌های شبه پارامتریک معمولاً ترکیسی از 
روش‌های پارامتریک و ناپارامتریک می‌باشند (۳۹ و ۲۱). 

روش مد نظر مطالعه حاضر روش ناپارامتریک شبیه‌سازی 
مونت کارلو است به این دلیل که از توزیع تجربی متغیر عملکرد 
براساس تغییرات آب و هوا در دوره گذشته برای اندازه‌گیری ریسک 


2 2106 ۷ ۷۷۳۵۵۲۵۲ -85 
16 2 -9 
دمک تمطاع ۳۱1 -10 
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استفاده می‌شود. مزیت عمده این روش آن است که هیچ فرضی را 
برای توزیع متفیر تحمیل نمی‌کند و اجازه می‌دهد براساس داده‌ها 
توزیع بدست آید اما قبل از آن باید شرایطی در نظر گرفته شود که به 
قرار زیر است: 

برای ارزیابی ریسک آب و هوا برای هر شاخصی در مرحله اول 
باید وابستگی بین شاخص مورد بررسی و ارتباط آن به آب و هوا 
تعریف شود. از آن جایی که در ۷۵ ریسک نامطلوب اندازه‌گیری 
می‌شود میزان ریسک در شرایط نامطلوب آب و هوا که ناشی از تغییر 
اقلیم برای دوره آینده است مد نظر قرار گرفته که ممکن است برای 
هر شاخص دیگری شرایط نامطلوب به نحوی دیگر تعریف شود. 

در این مطالعه ابتدا با استفاده از متغیرهای آب و هوایی در دوره 
پایه (۱۹۸۶-۲۰۱۵) ارتباط بین متغیرهای آب و هوایی و عملکرد گیاه 
رل کرک دعر رنه مس ود سرا تور آنن 
(۲۰۱۶-۲۰۴۵) نیز داده‌های آب و هوایی با استفاده از سناریوهای 
تغییر اقلیم حاصل از مدل ۲1200۷3 تحت سناریو انتشار ۸2 در 
مدل ریز مقیاس نمایی 1279-72 بدست می‌آید و در نهاییت مانند 
مرحله دوره پایه (۱۹۸۶-۲۰۱۵) با استفاده از مدل آکواکراپ عملکرد 
گیاه با در نظر گرفتن اطلاعات واسنجی شده در دوره‌پایه شبیه‌سازی 
ی سود 

پس از تعیین وابستگی باید توزیع شاخص براساس وابستگی آن 
به آب و هوا که از مدل آکواکراپ بدست می‌آید تعیین شود. هر متغیر 
تصادفی دارای توزیع احتمالی است که با استفاده از داده‌های تاریخی 
متفیر مورد نظر می‌توان نزدیکترین توزیع را با استفاده از آزمون‌های 
آماری مربوطه انتخاب کرد. برای این منظور با استفاده از آماره‌های 
«کی- دو۱» «کولموگروف- اسمیرنف"» و «اندرسون- دارلینگ "» 
اقدام به بررسی انطباق توزیع تجربی متفیر عملکرد بر توزیع‌های 
نظری مشهور می‌شود تا نزدیکترین توزیع احتمالاتی هریک از 
متیزض ین مود( 

پس از تعیین توزیع با استفاده از روش شبیه‌سازی مونتکارلو 
اقدام به شبیه‌سازی ۱۰۰۰۰ مشاهده برای عملکرد می‌شود تا توزیع 
مورد نظر بر مشاهدات شبیه‌سازی شده منطبق گردد. پس از آن با 
محاسبه معکوس تابع توزیع مورد نظر در سطح احتمال معین مقدار 
ارزش در معرض خطر بدست می‌آید. 

شبیه‌سازی مونت‌کارلو: شبیه‌سازی مونت‌کارلو یک الگوریتم 
محاسباتی است و در دسته‌ای از الگوریتم‌های محاسبه‌گر قرار دارد که 
برای محاسبه نتایج خود بر نمونه گیری تکرار شونده تصادفی اتکا 
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می‌کند: به دلیل اتکای این روش بر محاسبات تکراری و اعداد 
تصادفی. روش مونت کارلو اغلب به گونه‌ای تنظیم می‌شود که توسط 
رایانه و بر پایه توزیع احتمالی» اجرا می‌شود. با استفاده از توزیعات 
احتمالاتی» متغیرها می‌توانند احتمالات مختلفی از نتایج مختلفی که 
رخ می‌دهد را داشته باشند. در طی یک شبیه‌سازی مونت کارلو 
ارزش‌های تصادفی از توزیع احتمالاتی متفیر تصادفی با یک پروسه‌ای 
که فرآیند نمونه گیری نامیده می‌شود بدست می‌آیند» هر مجموعه 
نمونه گیری شده یک تکرار نامیده می‌شود و خروجی حاصل از آن 
نمونه ثبت می‌شود که این کار صدها و هزاران بار تکرار می‌شود و 
نتیجه یک توزیع احتمالاتی از خروجی‌های ممکن است به طوری که 
نمونه‌های ثبت شده با توزیع احتمالاتی متغیر تصادفی (توزیع ورودی) 
به طور کامل انطباق یابند. از این رو توزیعات احتمالاتی یک روش 
واقعی تری از توصیف نااطمینانی در متغیرهای تجزیه و تحلیل ریسک 
هستند (۱۸). عامل مهم دیگر در شبیه‌سازی مونتکارلو روش نمونه 
گیری است. هر تکنیکی که به تکرار و مدت زمان بیشتری برای 
انطباق به توزیع ورودی نیاز داشته باشد کارایی کمتری دارد که در 
این مطالعه از تکنیک نمونه گیری مکعب لاتين" استفاده شده است. 
تکنیک مکعب لاتین نمونه‌های طراحی شده را به مطابقت نزدیک‌تر 
با توزیع ورودی مجبور می‌کند و در نتیجه برآوردهای آماری سریعتر 
(با تکرار کمتر از نمونه‌گیری مونت‌کارلو) به توزیع مورد نظر همگرا 
می‌شود. نکته مهم در نمونه‌گیری مکعب لاتین طبقه‌بندی توزیع 
احتمالاتی است که در آن منحنی تجمعی توزیع در فواصل برابر در 
فا احتبال ی زا باه شاد بکزارها تسش 
می‌شود سپس یک نمونه بدون جایگزینی به صورت تصادفی از هر 
فاصله یا طبقه توزیع مورد نظر انتخاب می‌شود. این روش به دلیل 
اجرای سریع‌تر (با توجه به تکرار کمتر) و نمونه گیری از هر قسمت 
کاراتر از روش نمونه‌گیری مونت‌کارلو است (۲۶). 

ارزش در معرض خطر با وجود مقبولیتی که در میان فعالان 
ریسک پیدا کرده دارای معایبی است که مهم‌ترین آنها به صورت زیر 
است: ۱- به دلیل عدم برخورداری از ویژگی انسجام* یک سنجه 
تمام عیار نیست ۲- با یک سطح احتمال مشخص فقط در مورد 
بیشترین زیان ممکن بحث می‌کند و زیان‌های بالاتر از خود را نشان 
نمی‌دهد (۲۲). 
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شکل ۱- ارزش در معرض خطر و ریزش مورد انتظار با فرض نرمال بودن توزیع 
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بر این اساس می‌توان از مدل ارزش در معرض خطر شرطی 
(۷۵8) يا ریزش مورد انتظار (۳5)" به عنوان میانگینی از 
زیان‌های محتمل استفاده کرد. این روش به خاطر اینکه می‌تواند 
دنباله پهن بودن توزیع را در نظر بگیرد نسبت به ۷۵5 معیار بهتری 
برای ریسک است که در رابطه (۷) نشان داده شده است" 


۳ 1 ون ]وم - 0۷ 
0 0۳ 


در روابط بالا "۳ معکوس تابع ۳ (226 ارزش متفیر 
تصادفی در میانگین شرایط آب و هوایی و 0 سطح احتمال می‌باشد. 

در شکل (۱) ۷۵5 و 0۷15 با فرض نرمال بودن توزیع متفیر 
تصادفی در هر دو حالت نشان داده شده است. همانطور که شکل 
نشان می‌دهد قسمت زیان در سمت چپ توزیع مورد نظر است. 

آمار و اطلاعات مورد نیاز این مطالعه مربوط به متفیرهای 
هواشناسی (بارش, دمای حداقل و حداکثر, تابش خورشیدی» رطوبت 


اون 2 مباله ۷ لقطمتانلجمن -1 

2- ۳06۵۵۵08 11 

۳- در هر دو روش ۷216 و 216 ۷/) انحراف از حالت میانگین به عنوان معیار 

ریسک به دلیل محاسبه اختلاف عملکرد در شرایط نامساعد آب و هوایی از عملکرد 
در شرایط میانگین آب و هوایی محاسبه شده است. 


نسبی حداقل و حداکثر و سرعت باد) و بصورت روزانه می‌باشد که در 
دوره ۱۹۸۶-۲۰۱۵ توسط اداره کل هواشناسی استان آذربایجان شرقی 
)٩(‏ گردآوری شده است. همچنین اطلاعات مربوط به عملکرد گندم 
دیم از دوره ۱۳۷۰-۹۴ از مدیریت جهاد کشاورزی شهرستان اهر (۱۴) 
بدست آمد. 


نتایج مطالعه نیز مانند روش تحقیق در سه بخش الف) تولید 
سناریو اقلیمی» ب) شبیه‌سازی عملکرد و ج) اندازه‌گیری ریسک ارائه 
شده است. 

الف) تولید سناریو اقلیمی: همان‌طور که گفته شد از مدل 
1275-72 برای تولید اطلاعات روزانه دمای حداقل و حداکثر بارش و 
تابش خورشیدی استفاده می‌شود که برای این منظور نرم‌اف زار 1275 
5 »۷ به کار گرفته شد. آزمون‌های آماری به کار گرفته شده برای 
مقایسه توزیعات احتمالاتی و ارزیابی قابلیت مدل شامل آزمون «», 
آزمون «کولم وگروف- اسمیرنف دو نمونه‌ای» (5-) «مجئور 
میانگین مربعات خطای استاندارد» ( 4۸/57 «مجئور میانگین 
مربعات خطای استاندارد نرمال شده» ( ۸۷۲/5۳ ) و «ضریب 
همبستگی» ( ) می‌باشد که نتایج حاصل از آنها در جداول (۱) و (۲) 
ارائه شده است. نتایج برای متغیر بارش به صورت تجمعی ماهانه و 
برای سایر متغیرها به صورت متوسط ماهانه گزارش شده است. 
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جدول ۱- نتایج آزمون کلموگروف اسمیرنف (16-5) و (4) برای متغیرهای بارزش؛ دمای حداقل و حداکثر و تابش خورشیدی ماهانه مشاهداتی و تولید 
شده در سال‌های ۱۹۸۲-۲۰۱۵ در ایستگاه اهر 
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شکل ۲- مبانگین ماهانه بارش تجمعی. دمای حداقل و حداکثر و تابش خورشیدی مشاهداتی و تولید شده برای دوره ۱۹۸۲-۲۰۱۵ ایستگاه اهر 
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جدول ۱ نشان می‌دهد که براساس مقادیر آماره «» و آزمون 
سال‌های ۱۹۸۶-۲۰۱۵ برای هر ماه بدست آمده در هیچ یک از ماه‌ها 
اختلاف معنی‌داری (در سطح احتمال ۵ درصد) در مورد میانگین و 
توزیع مشاهداتی و تولید شده متغیرهای بارش دمای حداقل و حداکثر 
و تابش خورشیدی وجود ندارده که این موضوع در شکل ۲ به وضوح 
قابل مشاهده است. 

جدول ۲- مقدار ضریب همبستگی. آماره ۳۱۷6 و ۱۳۲ 
را براساس متوسط اطلاعات ماهانه در دوره سی ساله (۲۶۰ مشاهده) 
مشاهده که از میانگین اطلاعات ماهانه در دوره ۳۰ ساله محاسبه شده 
است» بدست می‌آید. دلیل این امر عدم وجود توزیع نرمال برای بارش 


طی سال‌های مختلف متفاوت باشد. مقادیر آماره 1818۳ برای هر 
سه متغیر بیانگر وجود اختلاف کمتر در تولید داده توسط مدل است. 
مقدار خطای نرمال‌شده مدل برای متغیر بارش در محدوده ۰-۱۰ و 
دمای حداقل و حداکثر در محدوده ۱۰-۲۰ قرار دارد که نشان‌دهنده 
به ترتیب کارآیی عالی و خوب مدل است. همچنین مقادیر ضریب 
همبستگی نشان دهنده ارتباط قوی بین مقادیر تولیداتی و مقادیر 
مشاهداتی در متغیرهای بارشء دمای حداکثر و دمای حداقل است. 
لازم به ذکر است که چون در این مطالعه اطلاعات تابش خورشیدی 
در دسترس نیست و از اطلاعات ساعات آفتابی برای شبیه‌سازی تابنش 
خورشیدی استفاده شده است. بر این اساس نمی‌توان ضریب 
همبستگی و آماره‌های ۳15۳8 و ۱158 را برای تابش 
خورشیدی محاسبه کرد. 


جدول ۲- ارزیابی مدل 1275-۷2 با استفاده از ضریب همبستگی و آماره‌های 13۷۲5۳ و 1۷۲5 
۵ ,۱۲۰۱۷۲۹۲ 20 ,۱۸۱۷۲۹۲ 20 اصمت‌لهی ممناهاه تنم عصمته ۳۵06 ۱۲۷۵ -عه۱ 0۲ «متافتاوهء 1۳6 -2 12016 


متغیرهای آب و هوایی 
معا منبه۷ عافهنل) آماره 
دمای حداکثر دمای حداقل بارش «0اه)نج(۳۳6۵ عناعناهاه 
۰ ۱ ۱۷۱۱20۱۱۱ مبانا)هو مرجم ات۷ 
234 1.01 084 ۹۳( 
1324 1870 327 ۱۹۵ 
096 0917 098 1 


مأخذ: یافته‌های تحقیق 


وعصن‌صن موه کمک بعمتتا0 


براساس نتایج حاصله از آزمون‌های مربوطه (جدول ۱ و ۲) مدل با 
دقت بالایی قادر به شبیه‌سازی متغیرهای بارش دمای حداکثر دمای 
حداقل و تابش خورشیدی می‌باشد از این رو می‌توان از اين مدل برای 
پیش‌بینی متفیرهای آب و هوایی برای دوره آتی استفاده کرد. در گام 
بعدی از مدل 1275-72 برای پیش‌بینی متغیرهای مربوطه برای دوره 
۲۰۱۶-۴۵ با در نظر گرفتن سناریو اقلیمی ۲2001۷3 تحت 
سناریو انتشار ۸2 استفاده گردید. برای این منظور از داده‌های آب و 
هوایی بدست آمده برای هر دوره تبخیر-تعرق و نیز از نرم افزار 1710 
2 بهره گرفته شد. نتایج حاصل, مقایسه داده‌های 
مشاهداتی دوره پایه و آتی در جدول ۳ براساس متوسط سالانه هر ماه 
نشان داده شده است. 

یافته‌های پژوهش بیانگر آن است که متفیرهای آب و هوایی 
حاصل از مدل ۳12001۷13 تحت سناریو انتشار ۸2 در دوره ۲۰۴۸۵- 
۶ در مقایسه با دوره ۱۹۸۶-۲۰۱۵ افزایشی است. به طوری که 
میانگین ماهانه دمای حداقل و حداکثر به ترتیب ۰/۵۶ و ۰/۶۲ درجه 
سانتی‌گراده بارش تجمعی ۹/۱۵ میلی‌متر تابش خورشیدی ۰/۰۶ 
مگاژول بر متر مربع در روز و تبخیر و تعرق مرجع ۱۸/۴۸ میلی‌متر در 


سال افزایش می‌یاد. نتیج حاصل با مطالعات خلیلی اقدم و همکاران 
(۱۶)» سیاری و همکاران (۳۶) و آبابایی و همکاران (۱) مطابقت دارد. 
روند افزایش بارش متغیر بوده و بیشتر در فصل زمستان و بهار روند 
فزایشی دارد که همزمان با آواخر فصل رشند گندم می‌باشد. بیشترین 
میزان افزایش بارش معادل ۲/۹۳ میلی‌متر در ماه اکتبر رخ داده که 
شاد با ماههای عقر و آیان مان کیت کشیه ج پاش ها 
حداقل و حداکثر در بیشتر ماه‌های سال افزایشی است و در ماه‌هایی 
که کاهش دما رخ داده مقدار تغییر بسیار ناچیز می‌باشد (کاهش دمای 
حداکثر فقط در ماه دسامبر و در مورد دمای حداقل ماه‌های ژانویه» 
اکتبر و نوامبر می‌باشد). بیشترین کاهش دما مربوط به دمای حداقل 
در ماه ژانویه به میزان ۰/۱۴ درجه سانتی‌گراد می‌باشد و بیشترین 
افزايش را در فصل بهار و تابستان زمان بلوغ و برداشت گندم تجربه 
می‌کند. تنش حرارتی پس از گرده‌افشانی در محدوده دمای ۲۵ تا ۳۱ 
درجه به طور معنی داری موجب کاهش عملکرد می‌گردد اما در 
شهرستان اهر افزایش دما تا ۲۱ درجه سانتی‌گراد به ندرت رخ می‌دهد 
بنابراین افزايش دما موجب تنش حرارتی در گندم شهرستان اهر 
نخواهد شد. تابش خورشیدی بیشترین افزايش را در فصل بهار در ماه 


۱۳۸ نشربه اقتصاد و توسعه کشاورزی. جلد ۰.۳۲ شماره ۲. تابستان ۱۳۹۷ 


می که معادل با ماه‌های اردیبهشت و خرداد می‌باشد تجربه می‌کند 
که این نیز برای رشد گیاه یک مزیت به حساب می‌آید به این دلیل 
کنیع کر ار ]رانا یعس با مت ان وم 
خورشیدی جذب شده دارد و میزان جذب تشعشع خورشیدی ارتباط 
مستقیم با میزان سطح برگ مزرعه دارد. اگر گیاه رشد پاییزه خوبی را 
گذرانده باشد و به اندازه‌ی مناسبی از رشد پوشسش گیاهی برسد در 
رشد بهاره می‌تواند تشعشع خورشیدی زیادی را جذب کند و به 
عملکرد خوبی دست یابد. بیشتر تغییرات تبخیر و تعرق در ماه‌های 
گرم سال (فصل بهار و تابستان) رخ داده که با نتایج حاصل از تغییر 
دما هم‌خوانی دارد. طبق رابطه (۲) تبخیر- تعرق با عملکرد محصول 
رابطه معکوس دارد یعنی با افزایش آن عملکرد روند نزولی خواهد 
داشت. بیشترین میزان تغییرات تبخیر در ماه ژوئن (خرداد و تیر) زمان 
برداشت محصول رخ داده است و کمترین میزان تبخیر-تعرق در 
فصل پاییز (ماه‌های نوامبر و دسامبر) است که سبب می‌شود 
محدودیت محصول از نظر دسترسی به مقدار آب ناشی از بارندگی 
کاهش یابد. به طور کلی میزان افزايش دمای حداکنر, دمای حداقل» 
بارشء تابش خورشیدی و تبخیر و تعرق مرجع به ترتیب ۳/۷۸ 
۷ ۰۳/۱۵ ۰/۴۱ و ۱/۴۷ درصد متوسط دوره ۱۹۸۶-۲۰۱۵ 

ب) شبیه‌سازی عملکرد: پس از واسنجی عملکرد برای دوره 
پایه» عملکرد با استفاده از داده‌های روزانه متغیرهای آب و هوایی 
برای دوره پایه و آتی شبیه‌سازی شد. برای دوره پایه متوسط عملکرد 
۰ تن بدست آمد در صورتی که متوسط عملکرد واقعی ۰/۹۵۴ تن 
می‌باشد. دلیل این تفاوت آن است که آکواکراپ عملکرد را با توجه به 
اطلاعات آب و هوایی شبیه‌سازی می‌کند و شامل مشکلاتی از قبیل 
آفات» بیماری‌ها و علف‌های هرز نمی‌شود. این در صورتی است که 
زنکا تیم یکی ازاشکاات تولین گننم حطقه اهر مس بابرا 
دوره آتی نیز عملکرد ۰/۹۹۹ تن بدست آمد که دارای روند افزایشی 
است. دلیل اين امر افزايش غلظت کربن دی‌اکسید و بارش برای دوره 
آتی می‌باشد. با توجه به اینکه روند بارش تحت سناریو ۸2 برای 
منطقه اهر افزایشی است و گیاهان دیم بیشترین وابستگی را به میزان 
بارش دارند افزایش عملکرد تحت این شرایط در نتیجه تغییر اقلیم 
امری طبیعی است. البته علاوه بر بارش افزایش غلظت کربن دی 
اکسید با افزایش عمل فتوسنتز به رشد گندم که در واقع جزو گیاهان 
4 است کمک می‌کند که این به عوامل دیگری همچون مقدار آب 
و نیتروژن بستگی دارد. 

ج) اندازهگیری ریسک: پس از شبیه‌سازی عملکرد برای دوره 
پایه و آتی ریسک عملکرد محاسبه گردید. بدین منظور ابتدا توزیع 
عملکرد بدست آمد سپس با استفاده از ۱۰۰۰۰ داده شبیه‌سازی شده 
به روش مونت کارلو مقدار ریسک با استفاده از روش‌های ارزش در 


معرض خطر و ارززش در معرض خطر شرطی بدست آمد. 


برای بدست آوزفن توزیع عملکرد از داده‌های شبیه‌سازی شده 
برای دوره پایه و آتی استفاده شده است. دلیل اين امر در نظر گرفتن 
«کی- دو» «کولموگروف- اسمیرنف» و «اندرسون-دارلینگ» اقدام 
به بررسی انطباق توزیع عملکرد با توزیعات نظری مشهور شد برای 
این کار نرم‌افزار 5.5 ۳16 ۳25۷ استفاده شد. دو نمونه از نزدیکترین 
توزیعات برای عملکرد توزیع بتا! و ویکبای " بدست آمد نتایج حاصل 
از آزمون‌های مذکور در سطح احتمال ۵ درصد در جدول ۴ نشان داده 
شده است. 

ای اه راشای تا انس نخان آسازفای 
مذکور نشان می‌دهد که براساس هر یک از آزمون‌ها فرض صفر مبنی 
بر هم‌ارز بودن توزیع مشاهداتی با یکی از توزیعات نظری مشهور که 
در اینجا توزیع بتا و ویکبای در نظر گرفته شده, رد نشده است یعنی 
ترش نا مرویگیای رسای ماس رای بسک ما شنی الاک 
۵ ِ ۴ 
نظیر بوکوشوا" (۴)گودوین" (۱۰) و پیش‌بهار و همکاران (۲۸) 
ویکبای را توزیع مناسبی برای عملکرد معرفی کرده‌اند در حالیکه 
نلسون * (۲۵) و چن و میرندا" (ع) در تحقیقات خود برای عملکرد از 
توزیع بتا استفاده کرده‌اند. 

پس از تعیین توزیع ۱۰۰۰۰ مشاهده برای هر یک از عملکردهای 
دو دوره با استفاده از توزیعات بتا و ویکبای شبیه‌سازی شد و ارزش در 
معرض خطر و ارزش در معرض خطر شرطی به صورت انحراف آن از 
میانگین با استفاده از نرم‌افزار 517ن2(8) محاسبه شد که نتایج حاصل 
در سه سطح احتمال ۱ ۵ و ۱۰ درصد در جدول ۵ نشان داده شده 


است. 


(- شعیر تصادقی: 2 فازای توزيع شا ببا بارامترهای یب استت: ار ایغ 
چگ الی اختمال آن بذصورت ۰ 23۳ 2۳ . باس درتالن کسه 
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اندازه گیری ریسک عملکرد ناشی از تغییر اقلیم در گندم دیم شهرستان اهر ... 


۱۳۹ 


جدول ۳- میزان و روند تغییرات میانگین ماهانه بارش, دمای حداقل و حداکثر, تابشس خورشیدی و تبخیر و تعرق پیش‌بینی شده برای دوره آتی 


(۲۰۱۳۲-۲۰۵) در مقایسه با دوره پایه (۱۹۸۲-۲۰۱۵) 


0۲ ,ر0۰۵)۱0۲6 66۳۵ را هه ره رصم رمااهازمنمع 0مامن0عن متاهزنه۲۱ «اطمرمصه مط ۵۶ ممهنه2۲ 20 1۳۵۵۵0 -3 12016 
(2011-2015) 06۵۲۱0۵0 »0261 6۳6 ۲۵ 0ع۲وصنصرمع (2016-2045) مهن مان مط) باه مت)ه ان مهم مرج 290 صمتاهن0 ۲2 


بارش 


21 3 ماه 
2 3 فٍ طاصما۷ 
عم 2۳ ۳ 7 ۱۳ و 
3 ۲ دسامپر نوامبر اکتبر سپتامبر آگوست جولای ژوئن 
۹۹ ۲( 0 م5 ناش انا[ صبال 
915 033- 37 393 246 049 5- 016 
افزایشی کاهشی کاهشی افزایشی افزایشی کاهشی کاهشی افزایشی 
ترتب ۱ ات۱ از ۱۱ ات۱۱۱۵ تومع1۱0 ات۱۱ ات۱ ات۱۱ 
1 1 4 4 1 4 1 4 
056 001 001- 4- 022 033 097 116 
زا فزایشی کاهشی کاهشی فزایشی افزاث فزایشی فا 
عممعموز حتعمععع1 حتعدع:0ع1 ۲۵۵6۵5 ۱۱۹ ات ۱۳ متعمععع1 عتعمعع 1 
3 1۹ 4( ‌ 
063 8.- 016 003 013 026 118 130 
زا کاه فا زا فزایشی زا فزایشی فا 
ممعموز متعمعتعع1 صتعفععم1 ات ۱۱ عتعمععع1 ان ۱۱ ۱۱ صتعمعع 19 
ِ 2 3 3 ظ 3 ظ ظ 
06 011 4- 6- 1- 03- 0.07- 046 
افزای افزایشی کاهشی کاهشی کاهش کاهشی کاهش افزای 
عممعموز حتعمععع1 متعدع0۲ع1 ۲۵۵۲6۵5 تعدع0۲ع1 متعحع:0ع1 ص۵ع۵۳ع۲ حتعمععع1 
ِ ظ 3 3 04 3 ظ ظ 
1948 2 44 005 1917 336 35 1-14 
زا کاه کاهشی زا فزایشی زا فزایشی فزای 
عمهعموز ص۵ع۵۳ع۲ متعدع:هع1 ات ۱۱۳ ۱۱ ات ۱۱۳ متعمععع1 ت۱۵ 
ِ ۹ 3 ِ ۹ ِ 3 


۱ متغیر 
می آوریل مارس خوریه ژانویه علاهتنه۷ 
«ه ۸۳ ۱/29 ۱3 130 
میزان تغییرات 9 
026 157- 332 ور 012 ۳ و 
20 ۲ 
ارت ۳ ارت ارت فه ای ب 
افزایشی کاهشی افزایشی افزایشی افزایشی روند تغیبرات لّ 
ات۱۱۵ توهع:0ع1 ات۱۱۱۵ ات۱۱۵ ات۱۱ هب ,5 
0 ۷2۳۱260 
4 4 4 4 4 
میزان تغییرات 
119 110 11 06 014 ۱۳9 و و 
امه 1 ۲ 
افزایشی افزایشی افزایشی افزایشی کاهشی روند تغییرات 2 
متعمع‌مم1 صتعفع 12 عتعمع‌مم1 متعمع‌مم1 آده۵عظ ۳9 طِ 
لمع ممتاهنته ۷ 
3 ِ 3 17 
میزان تغییرات 
113 133 126 082 0 ۳ و و 3 
امه 3 ۳ 
افزایش افزای افزای افزایش ۵ و ل 
یت یت و ره ثابت روند تغییرات و 3 9 
حاعمععع] ۱۵ حاعمععع] حاعمععع1 ت 
اصفاعممن لمع ممتاهنته ۷ 
ِ 3 ِ ِِ 
میزان تغفییرات و 
08 042 028 003- 04 تاه 3 3 
امه 3 4 
افزایشی افزایشی افزایشی کاهشی افزایشی روند تغییرات 5 3 
صنعدع‌هم1 صنکدع‌معمز طنعمع عم مومع صنعهع‌مم1 زار 
لجع ممتاهنته ۷ 9 
ِ 3 ِ 0 3 
میزان تغییرات و 
341 127 035 044 027 مهن 8 7 
اصم‌جند: ِ ۲ 
افزاي افزایشی افزایشی افزایشی افزایشی روند تغییرات 2 ِ 
متعمع‌مم1 صتعمع »12 صنعمع‌هم1 صنعمع‌هم1 صتعدع عم 0 4 
لمع ممتاهننه ۷ حِ 
ِ 3 ِ ِِ 3 اک 


مأخذ: یافته‌های تحقیة 


وعصنل‌صت مه کمک بعمتتا0 


طبق جدول ۵ مقدار ریسک در هر دو مورد ۷۵3 و 0۷21 در 
منطو احتمالاکی لت کاهش هن بابک که این ستغله طییی اییخ 
چون هر چه سطوح احتمالاتی افزایش می‌یابد کوانتیل تابع به سمت 
میانگین توزیع حرکت می‌کند و اختلاف آن از ميانگین کمتر می‌شود. 
مقادیر 0۷26 از ۷26 در سطوح احتمالاتی خاص بیشتر است چون 
5 میانگینی از زیانهای محتمل است در صورتی که ۷2۳ 
بیشترین زبان ممکن را در نظر می‌گیرد که سبب می‌شود حالت 
انحراف ازمیانگین ۷۵1۴ بیشتر شود. در مقایسه دوره‌های آتی و 
پایه زمانی در سطوح احتمالاتی خاص مقادیر سنجه‌های مورد نظر 
در دوره آتی کمتر از دوره پایه است به این دلیل که میزان عملکرد 
در دوره آتی نسبت به دوره پایه افزايش يافته است. در واقع زمانی 
که ارزش در معرض خطر آب و هوا براساس مطالعه پریتین سلیر و 
همکاران" (۳۰) حالت انحراف از ميانگین ارزش در معرض خطر در 
نظر گرفته شود در دوره آتی در مقایسه با دوره پایه ریسک عملکرد 
ناشی از آب و هوا کاهش می‌یابد. این موضوع از تغییر متفیرهای 


1- 


بارش دمای حداقل و حداکثر تبخیر-تعرق و غلظت دی 
اکسیدکرین ناشی می‌شود. در مورد توزیع بتا مقلدیر تغییر در سطوح 
احتمالاتی ۱ ۵ و ۱۰ به ترتیب برابر ۰/۱۳۹ ۰/۸۲ ۰/۶۱ تن در 
هکتار (که معادل ۰۱۴۳/۵۷ ۸/۷ و ۶/۳۹ درصد متوسط عملکرد واقعی 
در دوره پایه است) می‌باشد, این در صورتی است که برای توزیع 
ویکبای کاهش ریسک در دوره آتی معادل ۰/۲۹۷ ۰/۸۱ و ۰/۲۱ تن 
در هکتار (که برابر با ۳۱/۱۳ ۸/۴۹ و ۳/۲۴ درصد متوسط عملک رد 
واقعی در دوره پایه است) است. در توزیع بتا کاهش می‌پابد. براساس 
سنجه ارزش در معرض خطر شرطی نیز مقدار کاهش ریسک در 
دوره آتی برای توزیع بتا در سطوح احتمالاتی ۰۱ ۵ و ۱۰ درصد به 
ترتیب برابر ۰۰/۱۵۸ ۰/۱۱۵ و ۰/۹۳ تن در هکتار (معادل ۱۶/۵۶ 
۵ و ۹/۷۴ درصد متوسط عملکرد واقعی در دوره پایه) بوده و 
این ارقام برای توزیع ویکبای معادل ۰/۴۰۳ ۰/۱۴۸ و ۰/۷۷ (معادل 
۴ ۱۵/۵۱ و ۸/۰۷ درصد عملکرد واقعی است) می‌باشد. نتایج 
حاکی از آن است که میزان ریسک بدست آمده وابسته به توزیع 
احتمالاتی مورد نظر برای عملکرد است که می‌تواند براساس توزیع 
در نظر گرفته شده ارقام متفاوتی را شامل شود. 
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جدول - بررسی انطباق توزیعات نظری با توزیع تجربی عملکرد برای دو دوره پایه و آتی 
سح فصحتعدط ۲0۲ جمناتطتطاعنل لم‌تشوری 0اعزنز ما ات مممتاهطتاعنت آی‌ناهمع 0 2000068 0 ۱۳۲682210 رد ماطد۲ 
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وا 
مأخذ: یافته‌های تحقیة 
وعصنل‌صت موه ععک بعمتتا0 
نتیجه گیری و پیشنهادها است به طوری که در اين مطالعه روند بارش و دماهای حداقل و 


حداکثر تابش خورشیدی و تبخیر و تعرق افزایشی بدست آمد که 


یک از متغیرهای آب و هوایی مد نظر در این مطالعه تا ۱ ۱ 
هر یک از متغیرهای آب و هوایی مد نظر در این نتاس هام‌هایتضتلفت هیزان وت تافر گذاری آن سخاوت ان 


ز رند که جهت و میزان تاثیرگناری آنها متفاوت ۱ 
ی بیشترین میزان تغییر بارش در ماه مهر وآبان همزمان با فصل کاشت 


اندازه‌گیری ربسک عملکرد ناشی از تغییر اقلیم در گندم دیم شهرستان اهر ... ۱۵۱ 


گندم صورت گرفت. بارندگی‌های پاییزه بیشترین تأثیر را روی عملکرد 
دانه گندم دیم دارند زیرا قوی‌ترین ارتباط عملکرد دانه گندم با 
راکهاش ماه مر بان امتت کنه مرها تجوانهزتی: را امل 
می‌شود» هر چقدر میزان بارش بیشتر باشد بذر زودتر جوانه می‌زند و 
باعث می‌شود محصول مراحل رشد خود را به خوبی طی کند تا بتواند 
به عملکرد قابل قبولی دست یابد. همچنین ذغیره و نباشت مواد 


غذایی در بافت‌های گیاهی در زمان رشد پاییزه اتفاق می افتد و به 
مرور افزایش می‌یابده در زمان حاکمیت سرماء باعث غلیظ شدن شیره 
گیاهی و نهایتاً مقاومت گیاه در مقابل انجماد و سرمازدگی می‌شود. 
توسعه و گسترش و رشد بافت‌های ريشه در مرحله رشد پاییزه شروع 
و توان جذب آب بیشتر شده و در نهایت مقاومت لازم را مقابل 
تنش‌های خشکی کسب می‌کند (۲۹) 


جدول ۵- سنجه‌های ریسک در سطوح احتمالاتی مختلف برای متغیر عملکرد گندم دیم اهر 
۵ ۱۲۱۱۵۵ ۵۱۱۲60 ترجه 0 وامع۱ و)تازصاودامص اصعممن0 صز می‌تاع]۱۷ »1651 -5 12016 


۴ به صورت انحراف از میانگین 
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۶ به صورت انحراف از میانگین 
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۱ درصد جح كِ 
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مأخذ: یافته‌های تحقیق 


وعصن‌صن موه کمک بعمتتا0 


دماهای حداقل و حداکثر در بیشتر ماه‌های سال افزایشی است و 
بیشترین تغییر را در فصل بهار و تابستان تجربه می‌کند» تابش 
خورشیدی هم بیشترین افزایش در فصل بهار تجربه می‌کند که 
افزایش هر سه مورد برای رشد گیاه مزیت به حساب می‌آیند به این 
دلیل که منطقه اهر جز مناطق کوهستانی و سردسیر می‌باشد افزایش 
دما موجب تنش حرارتی گیاه نشده بلکه رشد گیاه را افزایش می‌دهد. 
افزايش تشعشع خورشیدی نیز باعث افزایش سطح تغذیه گیاه شده و 
در رشد آن تاثیر مثبتی دارد. تنها عاملی که تأثیر منفی بر افزایش 
عملکرد دارد افزايش قدرت تبخیر جو ناشی از افزايش دمااست که 
سبب می‌شود میزان آب در دسترس گیاه کاهش یافته و روی عملکرد 


تأثیر منفی بگذارد. از طرفی دیگر افزایش غلظت گاز کربن دی اکسید 
نیز باعث افزایش فتوسنتز گیاه شده و روی عملکرد تأثیر مثبت دارد. 
بنابراین براساس آنچه که گفته شد عملکرد محصول در اثر تغییر اقلیم 
مدل ۴1200۷3 تحت سناریو انتشار ۸2 دارای روند افزایشی است 
که براساس آن ریسک عملکرد در اثر تغییر اقلیم کاهش يافت که با 
توجه به جهت تغییرات و میزان تغییرات هر یک از متفیرهای آب و 
هوایی که در بالا دکر شد قابل انتظار است. با توجه به مطالب مذکور 
چون یکی از مهم‌ترین عواملی که باعث افزایش ریسک عملکرد 
می‌شود افزايش تبخیر از سطح خاک است که برای کاهش تبخیر 
استفاده از خاک‌پوش‌های مصنوعی مثل پشم شيشه و لایه‌های 
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پلاستیکی و خاک‌پوش‌های طبیعی مانند لش و لاش‌برگ توصیه ‏ صورتی که گیاه رشد پاییزه را به خوبی طی کند ریسک کاهش خواهد 
می‌شود. همچنین در کاهش ریسک عملکرد تغییر اقلیم اقدامات یافت. از این رو به کشاورزان توصیه می‌شود با کاشت به موقع از 
مدیریتی مزرعه بسیار مهم است که یکی از اين اقدامات موّثر تنظیم کوتاه شدن دوره رشد و کاهش عملکرد ناشی از آن جلوگیری کنند. 
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0۴ ۱۷۲۵۵۱۵۲ 0۴ 062766 مطا 10۶ واطمجمده‌تنت0ع۲ مطا 

۶ 0۳۵1[ صمتصه .ععصفط ممصصتان لصا موی ۵۶ ومعتص1 متقصصتان .2007 0۰ تامحصصهک مه ,۷ تتفداط رم 1۵00610 
۵۵۰ 119 م۲ ۷1۸ صهلو:۳۵ 18) .3)1(:133-142 بطع27عوعع 02۵۵5 ۲1610 

۶ ات0[ صمتصدی .صقیا ما ممتام‌نلمي اممظ/ گنه هه معصقطاه عاقصصنات .2010 ,0۳۲ 11مصصهک 4ص رخ 16۵00610 
۵۹2(۰ 12) .8)3(:508-520 مط1۵86276 دمم ۲16۵۱0 

4 6 ,ونوه۲ ۱ 0۲0تصهاو بلج ۱۷ ۲0 فعطم۸۵0۵۲۵0۵ .2012 .2 تالا 2۳0 و2 ط11 و.۷ 200 .لا لا یک ۲2۴ وبا زر[ 
۵ ۲۲۵۵۱۵۵ اصمحاوع 1۳۷ 

6۵۰/۷۷۷۷۰۰۵ 26 ع1اهآته۷خ ومتافتاهه مسانه‌تتوه .2016 :12820 باتهم که واعتم۱۷ 

۶ ما0 صوصخ .ععصهتبافطا تعافحولل اتمه مه آهعن‌انه‌تنوه ما فصمتاج*مصع .2001 بامصعلع ۷ مصع ۱۷۲ ۱۷۲۱۲۵۵02 
3(:65-0۰) 83 رطع [وتالتاهتتوظ 
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اندازه‌گیری ریسک عملکرد ناشی از تغییر اقلیم در گندم دیم شهرستان اهر ... 


,۴۰ 6۲ل2ظ)۲۳۵0۵۵ آلعا ,رو ۷۵۷۵۲ ری مووعصیام؟ ضا ر.ل لهطاقک ر.۲ امصصصصت! بل تام ریگ طهال02طافصفوظ .نا 170ظ ارم 
نموه 10۲ وعلوعاه0د ممتادامم0ج فصه وم وامجمصصا متقصصناه مصتلاه۷]۵0 .2015 ٩.‏ طموه۷۷۵۵9۵-0۵11 20 رنظه 0000 
-543)۳۴(:1019 راصعصصماحص هام1 مطا ۵۶ ممممزمو لفنک)ه-2(۵۵ ۷ هه و0عمهمنامه مصئودا فنص 1 مه متصلته9. ظ1 
.1027 

۰ 1841012866 وم ۵ ممتاقانملق فطع صم فصمتام‌صتاويه لمعمتاتاهاتتاوتل ۵۶ ععصمیالص م1 ,1990 ۴۰ هن ممواع۲ 
71-۰ :(12)01 ,ومن۱۵0۵۵00 انمض ۵۶ لفصستاهز لهماممن ط)ز۱۱0 

که عاصهانهبد عمتمایصته مصه وتوبامصه عفن . فتاه . مصلونا. ما . تیان (2004) . بصمتجعممزم) . ۲۸118208 
.۰0۵19206 ۵۵۰/۱۱۱۷۱۷۷ 

نع ها ۱۷1۵1۵ ۱۷۲۵۲2۵ وم معصمط متقصصتاه ۵۶ ۲۶۲6۵5 .2015 ٩.‏ ههتندنتدنا مج ۱۷ طع0ه2صمصم‌تطهان ریط تقطه‌ططورظ 
-83 :(702 رطتهععوم۴ ومتصصمممعظ باتهم تقصیمژ .0242 مهم هه طمهمعممه متاممهمصممع لمتاد۵و ۵۶ ممتاهع011م2 
۳۵۲۹12۴۰ 18) ,106 

۲ صناتطم۲۲ (۵ظ۷۲) ممصهتنافص و۵ ۷۷۵۵0۵۵۲-۳286۵0 .2015 .هر 1201۲20 2 رت 1۳2501 .و تلعطظ ر.ظ تقطو‌ططوزط 
نموه صهتصویا تقصتمژ . بممتاجعتام موه طمدهمتممه قلناممن ‏ مصز ۱۳-۱۷ :واصامن ‏ طمصه/۱۷۲ طا )ق۱۷۷۲ 0عت۱۴۵ 
۳۵۲۹1۵0(۰ 18) .37-62 :(903 ,م90۵0 وع۵808عظ 

0 1 نموه تج مج ۱۷۷۵2۲ .2016 ویک م2120 لتلطاه مضه ریت طمله12نلهطک ویاط تناممطمو۴ 
اه نموه رصماهمنصهعنن0 متمانم‌زری۸ع-11204 موزنههانه مه ۷۷۵۵۲ مطا ۵۶ لقصمژ .ععمنعع: 10مهع گه فعمتاتصمم 
0 

0صح 065601106 ۵ طمهمعممه رمگتصنا ۸ تفت اه فلج تعطاهع ۷۷ .2015 .1 1۲0ظ ۱۱۵11 م2 ر.ل اتعطامک ر.ظ تملمطامه)اع:۲ 
10-8 :(ظ)343 باصهحصوما 7 طظ 101۵۲ معط ۵۶ ممممزم .ععصقطه متقصصئاه حصم فعافی؟ فحصممهاً ۹6600۵0۲۵1 همهم 

۵ )12اه ما امه تمعن تمتجتجمته ۴۸0 ,2011 .۲ ومتعرع۳ 20 رین).1 منوا ر.ظ ماتالعاگ رین م۴ 
۰ 20]6۲ 20 ,3.1 صمنوته ۷ رامیاصه۷ عععمتع1م ,۷/۵/۵۲ 1۵ ۲۵5۵0۴886 

2 ع20010 :عصاتممط موه ها )صمجصمممصمهه منامگاعم ۷اع۸ .2015 ٩.‏ تقطع0صههات) هه ر.ظ تاطو000مک مصصمصط۴ 
۵۹۵(۰ 12) .91-113 :6024 ,اصمصصمومصه]۱۷ متا0۲)۶0ظ 24 ومله»عمتوصط اممصمصز لمصمز اصتهتفممم فتاه ما۷ 
۶ ۸۲۱۵۵۴ ۷۷۵۹۲ ۵۶ معصقطه متقصصتاه معبانظ 0۶ )مدحرحطاً ۰1۳6 .2011 بخ عدعظ 20 ر.ظ عمنتتنام ر.ظ صفهااناگ ر.ظ علل‌تام 
1073-0۰ :(3) 21 ,عمط آم)مممصممز صط لهطام1 27۲۰و معتحتعانا کجععع] عمط و006 )۱۷۷۸ :71۵106 

-3)8(:117 ,ط0تعععمک ومتطمصمعظ لقصتاول ,0مصامجور متاهاتاصصته مامت ماصما( ما ۵۶ توبن مه ,2003 یه تاه 
۳۵۲۹12۴۰ 18) ,138 

کصعصصممهصفط لول طا کافتیم ملد لمممتاتصمن و۷ عافتتعمیا ۷ .2008 .5 ۲۲/۵9۵۷ هه ریت مصتحتمو ریگ صللقف رنه 
270-4 ,)۵562۲ فهمتاهعم0 ها وله۵11) ۲ .صمتاهمتصصتاوه 20 

۶ 0۵۵۵۵۲1508 ,2011 ۱۷۲۰ تصهطمک تصحوع۲ هه ره تصتهووما۳ ۲2۲۱0 رب آه۸ صومصصحظ ریم طع220تله ولا تتهرهو 
تعلص وا تمه وم هه فعمامصدنهم ماقصصتاه ۵ ممتامن0ع۳۳ فطا ر۶0 (000۷۲2 هه ۳۸00۱۷۲3) واع0مه 00۷ #۷۵ 
۳6۵۳91280 18) ,912-925 :(25)4 بامق 20 عماه/۷۷ 0۶ فتاه .رصتعهها 28000صمفک :تناو عفوع) معصقطه مافصتاه 
60۰ طوناع ۲ 

۲ معافصصتاه طا م۱۵ 10۲ تماحتمصهع تمطا۱۷6۵ عتاففطههاو ۸ ۱۷۷۵-صکلضصا .2002 ,۲.۱۷۲ ۷«متبوظ مه رم ۷معمعصم5 
1-۰ ,1۵562700 60)عحصصقط) م۲ .0 .103و ۲۷ رأهتاصهوه و عون روعنل‌ن0او 

۱ ۷۵۲۱۵016 عامصت ۵۶ مصجمد همتامن0ع۳ هه ط12)10نصرزگ 2014 ۷۲ تامطع02مطک هه ,بط ۷۵۲۵0 ورگ تصصمتتجصمو تمم 1 
۶ هیام .رام‌طمطووها ماه هه مناد عقوع) واع۱۷]۵0 ممت)ههعت تجطامات مهم ناو مصنا تاد متفه 
۳۵۲۹12۰ 18) ,12028(:1-15 باصمصصمزن ص هه صقصت ۲ 

,16516 ۵۵6 .واه مره :وام0مصه وتاتهاه۷ متافمطه‌ماد و۷ آن که دج متالة ۷ .2012 ۷۰ 169270۷2 
۰ هم ط080عظ ۵۶ ماباتافصا روه‌عممته و 9۵02 ۵۶ هام۲۵ رمنامه:ظ ها والفه نوا دماعف 
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